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Abstract 

We present a detailed examination of thermalization after inflation for perturbative inflaton 
decay. Different interactions among particles in the plasma of inflaton decay products are considered 
and it will be shown that 2 — ^ 2 scatterings and particle decay are the important ones. We show 
that thermalization occurs after decays dominate scatterings, and that depending on the typical 
mass scale of inflaton decay products, different situations may arise. In particular, thermalization 
may be delayed until late times, in which case the bounds from thermal gravitino production on 
super symmetric models of inflation are considerably relaxed. We will also consider the case where 
the observable sector consists only of the MSSM matter content, and point out that flat directions 
with large vevs may result in earlier thermalization of the plasma and push the reheat temperature 
towards its upper limit. 



1 Introduction 



According to inflationary models 0, which were first considered to address the flatness, isotropy, 
and (depending on the particle physics model of the early universe) monopole problems of the hot 
big-bang model, the universe has undergone several stages during its evolution. During inflation, the 
energy density of the universe is dominated by the potential energy of the inflaton and the universe 
experiences a period of superluminal expansion. After inflation, the coherent oscillations of the 
inflaton, which behave like non-relativistic matter, dominate the energy density of the universe. At 
some later time these coherent oscillations decay to the fields to which they are coupled, and their 
energy density is transferred to relativistic particles, the reheating stage, which results in a radiation- 
dominated FRW universe. After the inflaton decay products are thermalized, the dynamics of the 
universe will be that of the hot big-bang model. 

Until recently, reheating was treated as a perturbative, one particle decay of the inflaton with 

i 

the decay rate T d (depending on the microphysics) , leading to the simple estimate T R ~ (T d M P i) 2 
for the reheat temperature (assuming instant thermalization) |2j . Tr should be low enough so that 
the GUT symmetry is not restored and the original monopole problem is avoided. In supersym- 
metric models there are even stricter bounds on the reheat temperature. Gravitinos (the spin-| 
superpartners of gravitons) with a mass in the range of 100 GeV-1 TeV (in agreement with low 
energy supersymmetry) decay after the big-bang nucleosynthesis. They are also produced in a 
thermal bath, predominantly through 2 — > 2 scatterings of gauge fields and gauginos. This re- 
sults in the bound Tr < 10 8 — 10 10 GeV, in order to avoid gravitino overproduction which would 
destroy the successful predictions of big-bang nucleosynthesis || . It has recently been noted that 
the initial stage of inflaton decay might occur through a complicated and non-perturbative process 
called parametric resonance, leading to an out-of-equilibrium distribution of final state particles 
with energies much higher than the inflaton mass (the preheating stage) [f|, |!|. This may have 
important cosmological implications, including gravitino overproduction which is most relevnt to 
our discussion. Gravitinos can be produced either through the scattering of particles in the preheat 
distribution 0, or directly during inflaton oscillations [[7], ||. |K|. In particular, it has been recently 
noted that the second mechanism can lead to efficient production of helicity ±| gravitinos, which 
possibly dominate the abundance of thermally produced gravitinos by several orders of magnitude 
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§0,0 0,0. 

However, it is believed that in a wide range of realistic inflationary models the final stage of infla- 
ton decay occurs in the perturbative regime and that the reheat temperature is determined therein. 
Thus a correct treatment of this stage and the subsequent thermalization of decay products is nec- 
essary. After all, parametric resonance may not occur or may be modified, and a mechanism may 
be found for the dilution of relics which were produced in dangerous abundances during preheating. 
In any case, the thermal production of such relics, gravitinos in particular, must be in accordance 
with nucleosynthesis bounds. 

Our aim in this letter is to give a more detailed examination of thermalization after perturbative 
decay of the inflaton. In this regime, inflaton decay products have energies higher than the mean 
thermal energy (if they were instantly thermalized) and a number density lower than the thermal 
equilibrium value. This implies that thermalization occurs when the interactions which increase 
the number of particles are at equilibrium. We will note that the inflaton decay is not sudden, but 
rather a prolonged process which starts very shortly after the end of inflation. We then consider 
the energy distribution of particles in the plasma of inflaton decay products, and compare the rates 
for their decay and 2 — > iV particle scatterings (with N > 3), as well as 2 — > 2 particle scatterings 
which only redistribute the energy among the scattered particles. It will be shown that (cascade) 
decay is the dominant process which leads to thermalization. We find that, in general, a seed in 
the plasma of inflaton decay products is thermalized first and then the bulk of the plasma will be 
thermalized very rapidly by scattering off particles in the thermalized seed. It is emphasized that 
the thermalization time depends on the the typical mass scale of those inflaton decay products 
which have interactions of gauge strength. Finally, we will consider the case where the observable 
sector consists only of the MSSM matter content, and discuss the effect of flat direction vevs on 
thermalization. 



2 The Energy Spectrum of Inflaton Decay Products 

After inflation the inflaton starts its coherent oscillations around the minimum of the scalar poten- 
tial, and at some later time it will decay to the fields to which it is coupled. In the perturbative 
regime the decay occurs over many oscillations of the inflaton field. This means that the oscillating 
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inflaton field behaves like non-relativistic matter consisting of a condensate of zero-mode bosons 
with the mass rrij. The inflaton decay rate Td can then be calculated from the one particle decay 
channel of these bosons. The value of Td depends on the nature of final state particles and their 
couplings to the inflaton. For example, if the inflaton has a Yukawa coupling h to (practically) 
massless spin-| fermions Yd = ^ m i'i while for the inflaton with gravitationally suppressed coupling 



to matter Yd ~ ■ Efficient inflaton decay happens at td = 3^ when H ~ T^. At this time the 
bulk of the energy density in the coherent oscillations of the inflaton is transferred to relativistic 
particles with an energy of order mj. From then on the universe is radiation-dominated, and the 
energy of relativistic particles is redshifted as a -1 oc t~^; where a is the scale factor of the universe. 
After the thermalization of inflaton decay products is completed the familiar hot big-bang universe 
is restored. 

However, inflaton decay does not suddenly happen at td- It is rather a prolonged process which 
starts once the inflaton oscillations start at t — £/, when H ~ Hi (Hi is the Hubble constant at 
the end of inflation). The comoving number density of zero- mode inflaton quanta at time t, n c (t), 
obeys the relation n c (t) = ni exp(— T^t); where nj is the inflaton number density at tj. The decay 
products have their largest number density at the earliest times, but they constitute only a tiny 
fraction of the energy density of the universe for t < td- The inflaton decay becomes efficient at the 
time td, and at this time the bulk of the energy density is carried by the decay products. Recently it 
has been noted that consideration of those decay products which were produced before td can have 
important implications for thermalization gravitino production []10| |, Affleck-Dine baryogenesis 



14f , and electroweak baryogenesis |fL5 



For ti < t < td the universe is matter-dominated and hence H — ^. In the time interval between 
t and t + ^H^ 1 = 2t, inflaton decay products are produced and their physical number density n at 
the end of this interval is 

H 2 

n = m[exp (-Y d t) - exp (_2r rf *)](— ) . (1) 

Here we have used the fact that for a matter-dominated universe the redshift factor for the 
number density of particles is (^-) -3 = (^§7) 2 (^ ne Hubble constant at the end of interval is y); 
where a/ is the scale factor of the universe at the end of inflation. The redshift factor for the 
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momentum of a particle (equivalently the energy of a relativistic particle) is (^) = (jfij) 3 - F° r 
t < y the above expression becomes n ~ f^~^ ri J- The energy of these particles is between (|) 3 m/ 
(for the particles produced at the beginning of interval) and mj (for those produced at the end 
of interval), and can be practically taken as mi for all particles. These particles have the highest 
energy in the spectrum of inflaton decay products at time 2t. On the other hand, particles which 
were produced in the interval [ti, 2ti] have the lowest energy in the spectrum. At 2t, their energy 

2 

has been redshifted (from mi) to (-^-) 3 m/, while their number density has been redshifted from 

After changing 2t back to t, and hence H to 2H, we find that for t < td the plasma of inflaton 
decay products consists of particles with energy E and number density n 

(jf)^mi < E <mi 

(2) 

3 -gfini <n< *-§4m 

such that 

3 

E z3HT d 

n — { — ttUi 3 

mi AHi y ' 

Particles with energy E Max = "mi and number density n Ma x = jjf^ n i are produced until t d . At 
this time, the inflaton decay is effectively completed, and almost all of the energy density is carried 
by relativistic particles. This implies that particles with energy mi are produced in the interval 
[y , td] for (practically) the last time. For H = T d , (2) gives the spectrum of inflaton decay products 
at t d : 

C^fYmi < E <mi 

1 " " (4) 
g£n, <n< §fcn, 

From then on the universe is radiation-dominated, (practically) no more particles with energy 
mi are produced by inflaton decay, and the energy and the number density of all particles will be 

1 3 

redshifted as t~? and t~? respectively. The spectrum of inflaton decay products for H < T d is 
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(5) 

(g)lg£„,<„<(g)ig£„, 

and we also have 



H*, E J 3F 



2 



Let us also find the occupation number of particles, as a function of their energy E in the 
spectrum. It is seen from (2) and (5) that particles with energy E at time t were produced during 

3 

a short interval At, at the time when the Hubble constant was H p = ( n E 1 ) 2 - The particle momenta 
at the time of production were in the range mj — miH p At < p < mi, and their number density at 
that time, from (1), was An ~ njT d A t(j^-) 2 . This results in f E oc 22 ~f (recall that / ~ and 
here Ap = miH p At). It is clear that Je does not change in time since An and p 2 A p are both 
redshifted as a -3 . 

3 Interactions in the Plasma of Inflaton Decay Products 

It is clear from (3) and (6) that the number density and the energy density of the plasma is 
completely dominated by particles with the highest energy Emo.x in the spectrum, both before and 
after efficient decay of the inflaton. Particles with lower energy have a considerably smaller number 
density and hence carry a much smaller energy density. It has been pointed out however, that lower 



energy particles in the spectrum may have an important role in thermalization [T3| . 

Thermalization is a process during which the energy density p of a distribution of particles 
remains constant, while their number density n changes in such a way that the mean energy of 
particles reaches its equilibrium value T. For a distribution of relativistic particles which consists 
of riB bosonic degrees of freedom and rip fermionic degrees of freedom at thermal equilibrium, we 
have p = ^{ns + \ n F)T 4 and n = ^-(ub + jnp)T 3 . Therefore, the ratios ^ and ^ determine 
the deviation from thermal equilibrium. If these ratios are greater than 0(1), the number density 
of particles should decrease, and hence the mean energy increases in order to achieve thermal 
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equilibrium. If they are less than 0(1), the number density should increase and the mean energy will 
decrease. In both cases, interactions which change the number of particles should be at equilibrium. 

Under the assumption that particles in the plasma decay very rapidly after scattering, thermal 
equilibrium is achieved when particle scatterings are efficient. If only scattering of particles with 
energy Emqx (which also have the highest number density riMax) is considered, the thermalization 
rate is 



2 n Max 

T T ~ a 2 (7) 

i^M ax 

-2\ 



where a is the gauge fine structure constant (an 0(10 ) number). After substituting for E Max and 
n Max from (5), we have 



HJ T r 2 



r>~« 2 (-) -^f-^nj (8) 
T d Hi mi 2 

Thermalization occurs when — H and is substantially delayed in general, resulting in a low reheat 
temperature which is consistent with the nucleosynthesis bound on the abundance of thermally 



produced gravitinos [16 



The scattering rate of particles with energy E s and number density n s off each other is ~ a2 j^z- 
Prom (6) we have -Jy oc E~%, which implies that particles in a seed with energy E s <C EMax are 
scattered off each other at a much higher rate and, therefore the seed could be thermalized much 
earlier than the bulk of the plasma. This has led to the notion of catalyzed thermalization of 
the bulk by a thermalized seed which happens if the number density of particles in the seed 
substantially increases after thermalization, and if particles in the bulk are rapidly scattered off 
these low energy particles. 

We have to bear in mind that 2 — > 2 scatterings do not change the number of particles and only 
redistribute the energy among the scattered particles. In order to increase the number density of 
particles in the plasma, the rate for one particle decay and/or 2 — > iV scatterings (with N > 3) 
should be at equilibrium. For a better analysis of thermalization it is therefore necessary to identify 
the relevant interactions and compare their rates. Here we list the important interactions of a 
particle with mass m and energy E s (assuming E s 3> m, a brief note on this will come later): 
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1- 2 — > 2 scatterings off other particles in the plasma. It is seen from (3) and (6) that n oc E? and 
hence particles with less energy have also smaller number density and energy density. This implies 
that the energy density of a seed with energy E s can considerably change only by scattering off 
particles with energy E > E s . For this to happen, it is also necessary that energy is redistributed 
among the scattred particles. Therefore we estimate the rate for scattering of a particle with energy 
E s off particles with energy E > E s and number density n, such that the transferred energy is 
AE > E s . The cross-section for this process is a ~ -j^, [], and the resulting scattering rate will be 
of order a 2 -^^. From (3) and (6) we have Jr oc E^ , which implies that the highest scattering rate 
is off particles with energy Emox m he spectrum. 

Therefore the rate for 2^2 scatterings which increase the energy of the particle E s is T scaU ~ 
a 2 0. If T scatt > H, particle will acquire an energy much greater than E s . However, the energy 
of the scatterers remains almost unchanged because their number density dominates over the other 
particles in the plasma. It is seen from (2) and (5) that Ts <g u oc frs (ia) before (after) efficient 
inflaton decay. Therefore if other interactions have negligible rates (compared with scatterings), all 
particles will finally have energies of order Euax- 

2 

2- Decay to other particles through kinematically accessible channels with the rate I" "decay ~ a it~ 
(-|r is the time-dilation factor). Such decays increase the number of particles and trigger a chain of 
cascade decays when r 'decay > H. If m 2 < scatterings in (a) increase the energy of the particle 
before it can decay. These scatterings will also result in a plasma-induced mass-squared of order 
a ™ Max PI. Therefore m 2 is at least of order a Z Max , even if the mass of the particle in the absence 

: Such scatterings occur as a result of effective contact interactions, i-channel processes or s-channel processes. 
Cross-section for the desired process, i.e. scatterings with energy transfer AE > E s , is readily found to be a ~ 
for a contact interaction (e.g., the quartic coupling of bosons coming from the D-term part of the action in 
supersymmetric theories). If E 3> E s , the cross-section for the t-channel processes with energy transfer E s < AE <C 

2 

E is much greater than However, these peocesses correspond to long range forces and are effectively screened 

2 

in the plasma. Overall, a ~ ■#-= is a reasonable estimate when all processes are taken into account. 

2 Since plasma contains particles of all energies in the spectrum we must verify that the inverse scattering is not 
important. In the perturbative regime of inflaton decay Je < 1, and the difference fE s fE Max — /e s +Ab/% m -A£ 
determines the direction in which the 2^2 scatterings proceed. As we found earlier Je E~^, which implies that 
scattering dominates over inverse scattering. 

3 Actually a 2 ^ is integrtaed over the whole spectrum. It is seen from (3) and (6) that the contribution from 
particles with energy Euax dominates the integral. 
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of plasma effects m is very small. As long as m 2 < a we have r decay ^ F scaU , and 2^2 
scatterings dominate decays. 

3- 2 — > iV scatterings off particles with energy £mbi and number density riMax, which happen at the 
rate T 2 ^n ^ a 3 J lMax T? 0. These scatterings increase the number of particles too. However, their 
rate is hierarchically smaller than those of 2 — > 2 scatterings and decays, and can be neglected in 
the analysis below. 

4 Thermalization and Nucleosynthesis Bound on the Re- 
heat Temperature 

The above arguments establish that the competition between 2^2 scatterings and (cascade) 
decays determine the thermalization time and thus the reheat temperature. Depending on the 
value of mo and parameters of the model for inflation nj, mi, and Td, different situations may arise. 
Our aim here is to study these different possibilities in detail. 

1- Consider the case where m 2 < a ™ A * { ax until very late times. In this case a particle with energy 
E s is scattered off particles with energy E Max before it can decay. Such scatterings increase the 
energy of the particle by an amount much greater than E s . This increase lowers the particle decay 
rate (recall that r 'decay E^ 1 ) and renders the decay inefficient. 
The scatterings become important when 

a 2 nMax > H (9) 

E m in E}^ax 

At this time particles with energy E m i n are efficiently scattered and acquire much higher energies. 
Subsequently the rate for scattering of particles with successively higher energies comes to equilib- 
rium, and they acquire much higher energies as time increases. All particles in the plasma have an 
energy of order EMax when 

2 n Max ^ tt /i n\ 

a — 2 > H ( 10 ) 

i^Max 



4 The extra factor of a appears because at least one more vertex with gauge strength interaction is needed, and 
there are more w- phase space factors from extra particles in the final state. 
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Shortly after that decay (now with the rate r decay ~ a 2 ^ ax ^ ) becomes effective, and a chain of 
cascade decays will lead to thermalization of the plasma. This yields H th ~ a 2 " Ma ^ , which after 
replacing for Euax and rtMax from (5), becomes 



and the reheat temperature is 



#tfc ~ 771 11 



This holds if m 2 < a^r 1 ^ for at least H > H th , giving 

m 2 < ——§—^ni (13) 

Therefore if the typical mass scale of particles in the plasma of inflaton decay products satisfy 
(13), the reheat temperature is given by (12). 

2- Consider the case where m exceeds the bound in (14). In this case mo 2 catches up with a ^ a3: 
at an earlier time, and m 2 ~ m 2 subsequently. This happens for (again we use (5)) 

H eq m (14) 

At this time T scatt ~ a 2 is at equilibrium for 



tii m I 

Particles with an (initial) energy less than E s in the spectrum have already had efficient scatterings, 
and therefore have acquired much higher energies. For particles with energies greater than E s decay 

2 

dominates over scattering, and is efficient when T decay ~ al2 E~ — ^- 

At H ~ H eq particles with energy E s , given by (15), decay efficiently and a seed is self- 

thermalized after a chain of cascade decays (]. The energy and the number density of particles 

5 Since the energy of decay products is smaller than the energy of the original particles subsequent decays in 
the chain occur in an even faster rate (recall that the decay rate is proportional to E^ 1 ). On the other hand, 
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in the thermalized seed are T ~ p s i ~ (E s n s ) J and T 3 ~ p s £ ~ (E s n s ) z respectively, and 
T 3 < ra A f a:r . Two processes are now competing with each other: particles with an energy- 
just above E s in the spectrum scatter off low energy particles in the thermalized seed at rate a 2 -^, 
and particles in the thermalized seed scatter off particles with energy Euax at rate a 2 ^Em* • ^ ^he 
first process is more efficient, particles in the bulk loose energy through scatterings off the thermal- 
ized seed. This increases their decay rate, triggers cascade decay, and leads to their thermalization. 
The thermalized seed grows very rapidly in this way until all of the bulk is thermalized. This is 



called the catalyzed thermalization of the bulk by a thermalized seed p3fl . If the second process 
is more efficient, particles in the thermalized seed acquire energies of order Em ax, and therefore 
later scatterings off the seed cannot lower the energy of particles in the bulk. In this case, cascade 
decay of particles with energies greater than E s (hence their thermalization) starts after the Hubble 
expansion has sufficiently redshifted their energy. Therefore the thermalized seed grows slowly until 
the first process catches up with the second one. At this time catalyzed thermalization occurs and 
the bulk becomes thermalized very rapidly. 

For the first process to be efficient two conditions are necessary. First, its rate should be at 
equilibrium 

a2 i~f - H (16) 

and, second, it should dominate the second process 

2 ^ 2 n Max , , -n. 

a — — > a — — (17) 

E S T TE Max 

where T = (n s E s ) J []. After substituting for E Max , nuax-, E s , and n s from (5) and (6), and using 
the fact that the seed becomes thermalized after a chain of cascade decays (i.e. when a 2 ^ 2 - — H), 
(16) yields 

^decay > T ' scatt regardless of the energy of the particle. Therefore scatterings are not important throughout the 
chain. At the end of the chain the number density of particles has increased while their energy has decreased. This 
implies that particles scatter off each other very efficiently, and therefore this can be defined as the time when the 
seed is thermalized. 

6 It is seen that (17) is satisfied when T 3 > riMax ■ This implies that catalyzed thermalization can occur 

when T 3 < n Max , and justifies the relation n s <C < nMax m above. 
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H < 

and (17) gives 



™18p 5^, 447 



n 2 m 28 fx 8 

h<{^A-) (19) 

Thermalization of the bulk occurs at the smaller H from (19) and (20), and for H < H eq (the seed 
cannot be thermalized for H > H eq because scattering dominates over decay then) [J 

„,18p 5 447 rJZim 28 tr 8 

H th = min [H eq , — 1 — , — r ~ r ) 20 

mfHj 1 d m i nr 

This holds if m is small enough such that rriQ 2 < a ™ Max before efficient inflaton decay, which, after 

using (4), gives 

nj < m < — 2 — n i ( 21 ) 



Therefore if m satisfies (21), H t h is determined from (20) and Tr ~ (H^Mpi) 1 . 
3- Consider the case where mn exceeds the bound in (21). In this case mn 2 dominates over a „ 
before efficient inflaton decay (i.e., for H > r<j). In order to find Hth, the above steps are repeated 
but the values of Euax, nMax, E s , and n s are replaced from (2) and (3), instead of (5) and (6). This 
leads to 

H « ~ar^T (22) 

The conditions for catalyzed thermalization of the bulk give 

H < ( ar "2 T > (23) 



7 In reality catalyzed thermalization of the bulk occurs after a number of efficient scatterings off the thermalized 
seed. It would be safer then to have a factor of 10 on the right-hand side of (16) and (17). However, since Tr oc Hth 1 > 
remains within the same order of magnitude even if changes by a factor of 10. 
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instead of (18), and 



,a 7 #/m 14 * 



m 7 r d m 

instead of (19). Therefore in this case 



H < ( ~ ^ 2 \ ) (24) 



n 9 V ?m 2 rir 2 3 /V 7 ff , 4 m 14 5 

= ™" K - J 3 °4' ) , 1 (25) 

mf tii ttiiL d ni 

It is also necessary that H t h > 1^ for thermalization to occur before efficient decay of the inflaton 
(otherwise we are back to case 2 in the above) which implies that mo must be indeed large. In this 
case the plasma is thermal when inflaton decay is completed at H ~ and Tr ~ (r^Mpj) 5 . 

In summary, depending on the value of mo and inflationary model parameters nj, rrij, and T^, 
we have 

,a 4 r/Mp| 2 i . i 

h-4 4 n /) ^ T R & i v dM Pl y (26) 
ii/ m/ 4 

Big-bang nucleosynthesis yields the bound Tr < 10 s — 10 10 GeV, which constrains the parameters of 
the model. This ensures that gravitinos are not produced in dangerous abundances in the thermal 
bulk. 

Even before thermalization of the bulk gravitinos are produced in the non-thermal bulk and 
the thermal seed. One may wonder whether the requirement Tr < 10 8 — 10 10 GeV also keeps 
the abundance of these gravitinos at a safe level. The number density of gravitinos which are 
produced in the plasma is in general proportional to n 2 (n is the typical number density of particles 
in the plasma) 0, which is much smaller before thermalization of the bulk. On the other hand, 
thermalization of the bulk releases a huge entropy and dilutes everything, including gravitinos which 
were produced in the thermal seed. In the case that thermalization occurs before efficient inflaton 
decay (case 3 in the above) it has been shown that the bound Tr < 10 8 — 10 10 GeV guarantees 
successful nucleosynthesis, despite the fact that the instantaneous temperature of the plasma is 
much higher at earlier times |TI|. We can expect that the same conclusion also holds for cases 1 



8 One factor of n comes from the scattering rate and the other one is from the number density of particles in the 
plasma. 
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and 2 above and therefore the strongest constraint on the model parameters is, in general, due to 
gravitino production in the thermal bulk. 

Two clarifying comments are to be made here. We have actually studied the thermalization of 
particles with mass mo- After their thermalization, the lighter particles reach thermal equilibrium 
through gauge-strength scatterings of these particles (e.g., two squarks exchange a gaugino and 
scatter to two quarks). Therefore only the typical mass scale of the sector is important for its 
thermalization and not the individual particle masses. Finally, the above analysis was done for 
relativistic particles and its consistency requires that Tr > mo- If Tr < mo, particles with mass 
mo become non-relativistic before thermalization occurs. These particles will then decay to other 
(relativistic) particles whose subsequent (cascade) decay leads to thermalization and determines the 
reheat temperature. The same analysis as above can be done in this case, but rrij and mo should 
be replaced by m and the mass scale of secondary decay products respectively. 

5 The MSSM Example and the Role of Flat Directions 

Here we consider a model which consists of the Minimal Supersymmetric Standard Model (MSSM) 
sector, the inflaton sector, and the low-energy supersymmetry breaking sector, where the latter 
two interact with the MSSM sector through gravitationally suppressed couplings. To keep the 
matter content minimal, no right-handed neutrinos or gauge sectors at the intermediate scale are 
introduced, and therefore the MSSM sector is the only one with gauge strength couplings. We 
consider a chaotic inflation model with the potential \m 2 (\) 2 for the inflaton cj), where the inflaton 
decays to other particles via gravitationally suppressed couplings (and hence there is no stage of 
parametric resonance decay in this model). From the COBE data on the isotropy of the microwave 
background radiation, mi ~ 10 13 GeV |L7|]- At the end of inflation cj) ~ 10 19 GeV which yields 



Hl ~ ^ ^ 10" GeV 
M Pl 

m = m I( j) 2 ~ 10 51 GeV 3 

3 

T d ~ ~ 10 GeV 

M Pl 2 
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Also for the MSSM sector, mo — 10 2 GeV. It is easily seen that with these values, (13) is satisfied 
and the reheat temperature derived from (12) is Tr ~ 10 6 GeV, which is far below the limit set by 
the nucleosynthesis bound. 

In the early universe however, it is possible to have much larger mass scales in the MSSM sector. 
Flat directions in the scalar potential of the MSSM (denoted as (p here) can acquire very large vevs 
during inflation, both in the minimal models Jl8| and in the no-scale models fl9| of supergravity. At 
the end of inflation ip ~ (i7/ M n ~ 3 ) ™~ 2 ||18|| , where n (not to be confused with the number density) 



is the order of nonrenormalizable superpotential term which lifts the flat direction, and M is the 
scale of new physics which induces the nonrenormalizable terms. The non-zero energy density of the 
early universe strongly breaks supersymmetry and induces a negative mass-squared of order —H 2 
for the flat direction. So long as other supersymmetry breaking sources do not induce a positive 
mass-squared of the same order, tp ~ (HM n ~ 3 ) ™ -2 and the flat direction vev decreases slowly. When 
other sources become dominant, the flat direction mass-squared becomes positive, oscillations start, 
and from then on ip is Hubble redshifted more rapidly. 

A non-zero vev for the flat direction breaks a subgroup of the MSSM gauge group. Those 
MSSM fields which have F- and D-term couplings to the flat direction acquire masses of order g<p, 
where g is a gauge or Yukawa coupling. This implies that in the early universe the mass scale of 
those fermions and scalars which are coupled to the flat direction, as well as the gauge fields and 
gauginos of the broken subgroup, can be much larger than 10 2 GeV. Thermalization in the MSSM 
sector can then happen earlier (case 2 or 3 in the above), leading to a reheat temperature higher 
than 10 6 GeV. For example, if mo ~ 10 10 GeV (which can be easily induced by the flat direction 
vev), thermalization occurs before efficient inflaton decay and Tr ~ 10 10 GeV. On the other hand, 
the flat direction vevs, and thus the induced masses, are rapidly redshifted after the flat direction 
oscillations start. The question is whether they remain large enough for a sufficiently long time 
so that thermal equilibrium is achieved. This depends on the time when oscillations start and the 
flat direction vev at the onset of oscillations. In the standard approach oscillations start at 
H ~ 777.3 ~ 10 2 — 10 3 GeV when the low-energy supersymmetry breaking takes over the Hubble- 
induced one. Recently it has been shown that many flat directions can start their oscillations much 
earlier, due to plasma effects |TJ]]. Therefore the mass scale of the model at early times, and hence 
the thermalization dynamics, has a dependence on the nature of the flat direction and its initial 
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vev. This may suggest another role for the supersymmetric flat directions besides baryogenesis 
and dark matter candidates ||21|| : they may also assist thermalization to happen earlier and lead to 
a higher reheat temperature. 

In conclusion the reheat temperature Tr may be below the upper limit (T^Mp;) 5 by several 
orders of magnitude, in which case the constraints on the parameters of the model are considerably 
relaxed f\ However, flat directions with large vevs may result in earlier thermalization and push Tr 
towards its upper limit. 



6 Conclusion 

We have considered thermalization after perturbative decay of the inflaton. In this regime the 

number density of inflaton decay products is smaller than the thermal equilibrium value, while the 

mean energy of particles is larger than its thermal value. Particles with the highest energy Emo,x 

dominate the number density and the energy density of the plasma of inflaton decay products, 

and 2-^2 scatterings off these particles increase the energy of other particles. We compared the 

rate for 2^2 scatterings, particle decay, and 2 — > iV scatterings and found that the first two 

are the important ones. Thermalization occurs when decays dominate scatterings and, depending 

on the model parameters, different situations may arise. We showed that the thermalization time 

and the reheat temperature has a dependence on the typical mass scale of inflaton decay products. 

If this mass scale is sufficiently large, a seed in the plasma is self-thermalized and the scattering 

of particles in the plasma off the thermalized seed leads to catalyzed thermalization of the bulk. 

Otherwise, the bulk is self-thermalized at late times, when all particles in the plasma have energies 

of order Emux- As a result, the reheat temperature can vary in a wide range. The strongest 

constraint (from nucleosynthesis) on the model parameters is derived for a very large mass scale when 

thermalization occurs before efficient inflaton decay. In general, a small mass scale considerably 

relaxes this constraint. We also considered the MSSM case as an example. It was shown that 

the mass scale of the model can be much larger in the early universe due to particle couplings 

9 For example, in models of Z3-tcrm inflation or supersymmetric models for new inflation the inflaton mass can be 
be as large as 10 15 — 10 16 GeV p2| . In such models the upper limit on Tr can be as high as 10 12 GeV. If mo is small 
enough, the reheat temperature will be far below its upper limit and hence in agreement with the nucleosynthesis 
bound. 
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to the flat directions with very large vevs. If the flat direction oscillations do not start early, 
these vevs stay large enough for a sufficiently long time. In this case, thermalization may occur 
earlier, thus leading to a higher reheat temperature and resulting in a stronger constraint on the 
model parameters. Finally, the dynamics of thermalization seems to be very detailed even in the 
perturbative regime of inflaton decay. 
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